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We show that a two-dimensional topological superconductor (TSC) can be realized in a hybrid 
system with a conventional s-wave superconductor proximity-coupled to a quantum anomalous Hall 
(QAH) state from the Rashba and exchange effects in single layer graphene. With very low or even 
zero doping near the Dirac points, i.e., two inequivalent valleys, this TSC has a Chern number as 
large as four, which supports four Majorana edge modes. More importantly, we show that this TSC 
has a robust topologically nontrivial bulk excitation gap, which can be larger or even one order of 
magnitude larger than the proximity-induced superconducting gap. This unique property paves a 
way for the application of QAH insulators as seed materials to realize robust TSCs and Majorana 
modes. 


PACS numbers: 73.43.-f, 81.05.ue, 71.10.Pm, 74.45.+c 

I. INTRODUCTION 

Majorana modes can naturally exist in topological su¬ 
perconductors (TSCs)ji“— The intrinsic TSC has been 
predicted to exist in superconducting Sr 2 Ru 04 with p- 
wave paring statei^ii However, this has not yet been 
experimentally confirmed. Recently, many efforts have 
been devoted to design artificial TSCs4“— So far, most 
studies focus on the effective p-wave superconductors in 
hybrid systems with conventional s-wave superconduc¬ 
tors in proximity to strong topological insulators^ semi¬ 
conductors with strong spin-orbit coupling (SOC))^“— or 
ferromagnetic atom chains Some attention has also 
been paid to the conventional s-wave superconductors 
coupled to quantum anomalous Hall (QAH) insulators 
such as topological insulators with magnetic dopants.— 
Among all the above TSCs, multiple spatially overlap¬ 
ping Majorana modes, which greatly benefit the trans¬ 
port properties, can only coexist in one-dimensional (two- 
dimensional) TSCs belonging to Class 
(D^“— ) with integer topological invariant— In reality, 
the one-dimensional TSCs in Class BDI can easily reduce 
to the ones indexed by Class D with zero or one Majo¬ 
rana mode— As for the two-dimensional TSCs in 
Class D, the number of the Majorana modes or the Chern 
number is limited upto two— More Majorana modes 
or larger Chern numbers are limited by large chemical 
potential (i.e., very high doping) and an overall much 
smaller bulk excitation gap than the proximity-induced 
superconducting gap—>2^ 

In this work, we show that a two-dimensional TSC can 
be realized in a hybrid system with a conventional s-wave 
superconductor proximity-coupled to a QAH stated due 
to the Rashba SOCSi and exchange field in single layer 
graphene. Interestingly, with very low or even zero dop¬ 
ing near the Dirac points, i.e., two inequivalent valleys, 
the TSC from the QAH state has a Chern number reach¬ 
ing as large as four, hosting four Majorana edge modes. 


More importantly, these Majorana modes are protected 
by a bulk excitation gap, which can be larger or even one 
order of magnitude larger than the superconducting gap 
from the proximity effect. This is in strong contrast to 
the case of effective p-wave superconductors where the 
excitation gap is always smaller than the superconduct¬ 
ing gap— As the large topologically nontrivial gap has 
been shown to be probably most important for applica¬ 
tions in topological insulators^i^ topological crystalline 
insulators,— and QAH insulators,i22ri^ our finding, i.e., 
reporting a large bulk excitation gap in the TSC is cru¬ 
cial to the field of TSCs and Majorana modes. This paves 
a way to obtain robust TSCs and Majorana modes us¬ 
ing the QAH states. We also address the experimental 
feasibility of the TSC from the QAH state. 

This paper is organized as follows. In Sec. H, we 
present our model and lay out the tight-binding Hamil¬ 
tonian of single layer graphene. Then, we calculate the 
topological invariant in Sec. HI. We further presents the 
results on the phase diagram, Majorana edge states and 
bulk excitation gap in Sec. IV. Finally, we summarize 
and discuss in Sec. V. 


II. MODEL AND HAMILTONIAN 

The real-space tight-binding Hamiltonian of single 
layer graphene with the Rashba SOC, exchange field 
and proximity-induced s-wave superconductivity is given 

- M X + W X 

ia. ioi 
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where (*,j) represents the nearest-neighboring sites and 
Cia (c|q) annihilates (creates) an electron with spin a 
at site i. The first term stands for the nearest-neighbor 
hopping with t = 2.7 eV— being the hopping energy. The 
second term denotes the Rashba SOC with A, cr and 
representing the coupling strength, Pauli matrices for real 
spins and a unit vector from site j to site i, respectively. 

in the third term is the chemical potential. 14 (A) 
in the fourth (fifth) term corresponds to exchange field 
(superconducting gap from the proximity effect). 

To start, we transform the Hamiltonian of Eq. Ill) to 
the Bogoliubov-de Gennes (BdG) one in the momentum 
space. Specifically, 

^=^E‘^k^^BdG(k)<l>k (2) 

k 


where = (V'At(’^)> V'btW, V'A;(k), V'bj.W, V'A 4 ,(-k), 
V’B 4 ,(-k), --0At(-k), -■i/'Bt(-k)) with t/’|„(k) creating 
an electron with spin a and momentum k counted from 
the momentum T at sublattice i (i = A, B) and 


i?BdG(k) = (^ ^ 


A ^ 


(3) 


M. The detailed calculation is shown in Appendix |B] 
Obviously, our system is topologically trivial in the case 
of |14| < |A|. As for |T4| > |A|, we have ten critical 
chemical potentials in order, i.e., /ii _2 = 3t± 

M3.4 = t ± - A^, - A4 ^7,8 = 

—t ± \/V^ — A^, and /.tguo = —3t ± \/V^ — A^ by as¬ 
suming |t4|, |A| <C t, which divide the system into eleven 
topological regimes. 

These topological regimes are characterized by the 
Chern number Ci since i?BdG(k) belongs to Class D with 
integer topological invariant)^ Ci can be calculated by^ 

= ^ [ d^kUyik) (5) 

Jbz 

with the Berry curvature 

/xy(k) ~ /n)wL(k)[9fe,i?BdG(k)]n„(k) 

m,n 

X uUk)[dk,Hji^GmuM/[EM - A„(k)]2. (6) 

Here, Mm(k) is the m-th eigenvector of i?BdG(k) with 
the corresponding eigenvalue being Am(k); /^ = 1 (0) 
for occupied (empty) band. The Chern number of all 
topological regimes is given by 


i?e(k) represents tight-binding Hamiltonian without the 
s-wave superconductivity, which can be written as 


He{k) 


( V, /(k) 0 h,{k)\ 

r(k) f4 h*{k) 0 

0 h2ik) -V, /(k) 

\hiik) 0 r(k) -14 / 


(4) 


where /(k) = —1[(2 cos ^ cos + cos-^) — 

i(2cos^sin^ - sin^)], /ii(k) = -A[(cos^ + 

■y3sin^)sin^ -I- sin-^] - iA[-cos^ + 

cos ^^(cos ^ + v^sin ^)] and /i 2 (k) = A[(\/3sin ^ — 

cos !f) sin ^ - sin + ^[cos ^ - cos ^(cos !f - 

•\/3sin^)]. Note that the lattice constant is set to be 
unity in the calculation for simplicity. 


1 (-1), /r2 < M < ^1 
-3 (3), /r4 < /i < /ra 
4 (-4), /re < M < M5 
-3 (3), /rs < /i < /ry 
1 (-1), Mio < M 
, 0 , other regimes 


(7) 


when 14 > 0 (14 < 0). It is seen that jCil = 1 (3) near 
the momentum T (M) point, which is consistent with 
the number of the zero energy states in Ref. |2^ These 
Majorana modes require very large chemical potential (of 
the order of eV), clearly unachievable experimentally. It 
is noted that the study on the Majorana modes near the 
Dirac points is absent in Ref. In this work, with 
very low or even zero doping near the Dirac points, i.e., 
K (two inequivalent ones), we have a Chern number as 
large as four. 


III. TOPOLOGICAL INVARIANT 

Before investigating the topological properties of 
77BdG(k), we first identify the gap closing conditions. 
The gap closing of the BdG Hamiltonian iJBdG(k) is 
equivalent to the existence of bulk zero energy states due 
to particle-hole symmetry. The condition for bulk zero 
energy states is obtained by calculating det(i/BdG) = 0 . 
We find that the gap closes at the momenta T (single 
one), M (three inequivalent ones) and K (two inequiva¬ 
lent ones) points with the corresponding conditions given 
by (/t ± 3t)2 = ^ ^ y-i _ ^2 

= 14 ^ — A^, respectively. It is noted that -|- (—) stands 
for lower (higher) energy band at the momentum T or 


IV. RESULTS 
A. Phase diagram 

In the following, we focus on the investigation near 
the Dirac points. We first study the topological phase 
diagram as shown in Fig. UKa). The phase boundaries 
between the topological and nontopological superconduc¬ 
tors (NTSCs) are determined by the dashed curves, i.e., 
U? = /i^ + A^. To further distinguish the TSCs (i.e., 
> fP + A^), we suppress the s-wave superconductiv¬ 
ity. Without the s-wave superconductivity, we show the 
bulk energy spectrum of the low energy effective Hamil¬ 
tonian near the Dirac points (see Appendix in 
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Fig. urb). When the chemical potential lies in the gap 
(|/i| < Eq), eg., /iin, the system behaves as a QAH state 
with the Chern number |A^| = 2.— Note that Eq is the 
absolute value of the minimum (maximum) energy of the 
conduction (valence) band with the formula given in Ap¬ 
pendix This QAH state in proximity to an s-wave 
superconductor becomes a TSC with the Chern num¬ 
ber 2|A^| = 44^ (see regime I). When the chemical po¬ 
tential is tuned out of the gap below the upper limit 
1141, eg., flout, the system is in a metallic phase with 
two Fermi surfaces in each valley as shown in Fig. [TJb) 
{K' valley is not shown here). With the s-wave super¬ 
conductivity included, the effective paring near each of 
these four Fermi surfaces is equivalent to that of a p-wave 
superconductori^iiiiSl Each of these effective p-wave su¬ 
perconductors hosts a Majorana edge mode, which is in 
agreement with the Chern number near the Dirac points, 
i.e., jCil = 4. This effective p-wave superconductor from 
metal is labeled as regime IF Similarly, the NTSCs (i.e., 
< fi^ + A^) can also be divided into two regimes, i.e., 
regime III (from the QAH state) and regime IV (from 
metal). 


B. Majorana edge states 


As the effective p-wave superconductors (regime H) 
have been widely investigated in the literature,— we 
concentrate on the TSC from the QAH state (regime 
I). The Majorana edge states are studied in thick 
graphene ribbons. The numerical method is detailed 
in Appendix [C] We plot the energy spectrum of zigzag 
graphene ribbon near the K and K' points in Figs. HJa) 
and (b), respectively. We find that there exist four zero 
energy states in each valley. These eight states can be 
divided into two categories, i.e., four propagate along 
the same direction +x {—x) determined by the group 
velocity Vx = > 0 (< 0). Moreover, the four 

states in the same category are at the same edge, which 
is in agreement with the magnitude of the Chern num¬ 
ber, i.e., \Ci\ = 4. This indicates that these eight zero 
energy states are topologically protected Majorana edge 
states. Specifically, we choose two of them at the same 
edge with Ua, > 0 near the K point and show the real 
space probability amplitude of the one with smaller and 
larger \kx\ in Figs. [5](c) and (d), respectively. Note that 
we separate the A and B sublattices by the blue solid and 
red dashed curves. It is seen that the amplitudes of both 
A and B sublattices in two Majorana edge states show 
obvious decay and oscillation. However, the penetration 
lengths are different between these two Majorana edge 
states. 


V 

A 0 (A = 0) 

s 

\ 

\ 

V 

s 

X 

X 

X 

X 

(a) 

TSC 

(QAH) 

I 

✓ 

III / 

II ^ 

TSC 

(Metal)^*' 

y 

y 

y 

y 

IV 






/ m 

SC \ 

NTSC 


/ (QA 

H) \ 

(Metal) 







s 

•X 

✓ 



X 

✓ 



X 

y 



S 

y 



N 

y 



X 

y 



X 

y 






FIG. 1: (Color online) (a) Topological phase diagram in the 
(/r, 14) space with A V 0 or A = 0. The dashed curves, i.e., 
V? = A^ -|- fi^ are the phase boundaries between the TSC 
and NTSC whereas the dotted ones, i.e, fi^ = Eq stand for 
the phase boundaries between the QAH state and metal, (b) 
Bulk energy spectrum of near the K point with ky = 0 
and A = 0. 14 ( — 14) is the upper (lower) limit of the chemical 
potential in the topological nontrivial regime (VI? > fi^ + A^). 
fiin and flout stand for the chemical potential in and out of 
the gap, respectively. Vz = 6 meV and A = 4 meV. 


C. Bulk excitation gap 

1. Chemical potential dependence 

The above Majorana edge states are protected by a 
bulk excitation gap of the TSC from the QAH state. 
With different chemical potentials chosen in the gap of 
a QAH system, the bulk excitation gap as a function 
of the proximity-induced superconducting gap is plot¬ 
ted in Fig. 131 In the A = 0 limit, the system can be 
considered as two copies of QAH insulators as shown in 
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FIG. 2: (Color online) (a) and (b) represent the energy spec¬ 
trum of zigzag graphene ribbon with the Rashba SOC, ex¬ 
change field and proximity-induced s-wave superconductivity 
near the K and K' points, respectively, (c) ((d)) Real space 
probability amplitude \%l)\ across the width for the Majorana 
edge state with a smaller (larger) \kx\ near the K point at 
one edge (i.e., y — 0) with Ua, > 0 (only part of the ribbon is 
shown). A (B) refers to A (B) sublattice. The fluctuations 
of I?/)I at the positions far away from the edge are due to nu¬ 
merical error. Here, 14 = 6 meV, A = 4 meV, y = 0 and 
A = 2 meV. 


FIG. 3: (Color online) Bulk excitation gap Sgap of the TSC 
from the QAH state as a function of A. The solid curves with 
diamonds, crosses and plus signs correspond to the numeri¬ 
cal results at /r = 0, 1 meV and 2 meV, respectively. The 
analytical results at /r = 0, 1 meV and 2 meV are separately 
represented by the symbols of squares, upward triangles and 
downward triangles. Note that for the analytical results at 
/i 7 ^ 0, only two limits, i.e., A 0 and A ~ Ac are calcu¬ 
lated. In addition, the dotted (dashed) curve corresponds to 
Sgap = A (i?gap = lOA). Vz = 6 meV and A = 4 meV. 


Fig. mb) but with an energy shift of —fi (/r) for the par¬ 
ticle (hole) one. Then, the bulk excitation gap of our 
system is determined by these two QAH insulators, i.e., 
Egup = Eq — 1^1. This nonzero bulk excitation gap in 
the limit A = 0 strongly indicates that the bulk ex¬ 
citation gap can be much larger than A especially for 
small A. It is emphasized that the nonzero excitation 
gap in the A = 0 limit is totally different from the case 
of the effective p-wave superconductors where the exci¬ 
tation gap is exactly zero in the limit A = 0.^- At the 
critical point Ac = \/'V‘^ — the bulk excitation gap of 
our system becomes zero. In between, the bulk excitation 
gap shows a monotonic decrease with increasing A. We 
emphasize that during this process, Figap can be larger or 
even one order of magnitude larger than A by referring 
to Figap = A (dotted curve) and Figap = lOA (dashed 
curve). For example, i?gap = 4.02 meV (A = 0.3 meV) 
at p = 0; Figap = 3.22 meV (A = 0.3 meV) at p = 1 meV; 
Figap = 2.23 meV (A = 0.2 meV) at p, = 2 meV. This 
marked enlargement of the gap is in strong contrast to 
the effective p-wave superconductors where the bulk ex¬ 
citation gap is always smaller than the induced supercon¬ 
ducting gapj^“— This makes our proposal, i.e., the TSC 
from the QAH state, very promising for the realization 
of robust Majorana modes in experiments. 

To have a better understanding of the behavior of the 
bulk excitation gap of the TSC from the QAH state, 
we also perform an analytic derivation. Near the Dirac 


points, the BdG Hamiltonian iFBdG(k) in Eq. ( 3 ) can be 
expanded as a low energy effective one with iFe(k) [see 
Eq. ([ 3 ])] being replaced by Ff|®(k) [see Eq. (lAll) ]. The 
secular equation of the eigenvalue E is det[iFBdG(k) — 
EIsxs] = 0 where Fgxs is a unit matrix. After a careful 
calculation, we have 

[aj - + 4ai(A|. - - pI4) -I- 4 a 2 (p^ - A|.)]^ 

- 64Hfa3(A^ - + 8 [aip - 2 (p -f Vz) 

X (m" - A|e)][(p + Vz)(a^, - 4H>3) - 2paia2] = 0 

( 8 ) 

with ai = a 2 — as + 04 , 02 = as = 04 = 

A^ — + p^, Vf = 3t/2 and Xr = 3A/2. Note that we 

focus on the calculation near the K = ( 47 r/ 3 , 0 ) (r = 1) 
and set ky = 0 hy considering the isotropy of the low 
energy effective Hamiltonian. It is very difficult to obtain 
the eigenvalues by solving Eq. ([5]) directly. Instead of the 
eigenvalues, we are interested in the bulk excitation gap 
here. Differentiating Eq. ([5]) with respect to 02 and then 
employing the extreme value condition of the excitation 
gap (i.e., = 0), we have 

ai - 920 il - gias - go = 0 (9) 

where p2 = 8(02 + 04) -|- 2 ( 2 A^ — p^ -I- 2 V^), gi = —3(02 + 
04)^ + 4 (— 2 A|j — V^ p^)(q;2 + a4) -|- 4 a 2 (A|j -I- p^) — 
8 V^{Xj^ + g?) and po = (a2 + 04)^ — 2(p^ — 2A|j)(a2 + 
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aiY + 4(q!2 + ctA)[a2{-X\ - - A^)] + 

8a2(A^"' — A|j) + 8a4A|.(^^ — A|j). 

At ^ = 0, Eq. dH) can be simplified to (02 + — 

a 3 )( 4 a§ + qiuz + 92 ) =0 with 91 = — 8(02 + 04 ) — 
16(A|j + V^) and 92 = 4 (a 2 + <a 4 )^ + 16Afj(a2 + ai) — 
16A|ja2+32A|jE/. Since the equation 4a|+ 9103+92 = 0 
is inconsistent with the gap closing condition, we only 
have a 2 + 04 — Q !3 = oi = 0. With this condition 
together with Eq. ([8]), one obtains the bulk excitation 
gap Egap = £’ 0(1 — lAj/lE^I), which is linearly depen¬ 
dent on A and agrees very well with the numerical re¬ 
sults as shown in Fig. [31 Specially, for £gap > |A| 
(£gap > 10|A|), we have |A| < £o|V;|/(£o + = Ai 

(|A| < Eo\V^\/{Eo + 10|44|) = A 2 O.lEo). These con¬ 
ditions will guide the experiments to obtain robust TSCs 
and Majorana modes. As for the case of ^ 0, it is very 

difficult for us to obtain an exact analytic solution. Only 
the analytical results in two limits, i.e., |A| ^ Ac and 
A ^ 0, are given. In the |A| ^ Ac limit, we have £gap = 
(Ac - \A\)Al\^^^ - A||/v^K ^(Ag+Aj,)(/r4 + Aj,Ag ). In 
the limit A ^ 0, £gap = with 

+ KXr + - V^)Eo and 

W2 = - m^)(A| + E/)£o - 16 A?jE/(A|, - - 

+ X\V^) / {X^i^ + V^) by assuming 0 < ^ < Eq. The 
analytical results at /x 7 ^ 0 in both limits agree fairly well 
with the numerical ones as shown in Fig. 3. 

2. Rashba SOC strength and exchange field dependences 

We then turn to investigate the effects of the Rashba 
SOC and exchange field on the bulk excitation gap of 
the TSC from the QAH state. In Figs. lUJa) and (b), 
we plot the dependence of the bulk excitation gap on 
the proximity-induced superconducting gap at /x = 0 un¬ 
der different Rashba SOC strengths and exchange fields, 
respectively. It is seen that the bulk excitation gap 
increases with the increase of either the Rashba SOC 
strength or exchange field. This can be easily understood 
from Fgap = £ 0(1 ~ l^l/|4^z|) mentioned above where Eq 
(see Appendix [^ increases with increasing Rashba SOC 
strength and exchange field. 

V. SUMMARY AND DISCUSSION 

In summary, we have proposed that in the presence 
of proximity-induced s-wave superconductivity, the QAH 
state due to the Rashba SOC and exchange field in sin¬ 
gle layer graphene can become a two-dimensional TSC. 
With very low or even zero doping near the Dirac points, 
i.e., two inequivalent valleys, we show that this TSC, 
which exhibits a Chern number as large as four and hosts 
four Majorana edge modes, has a bulk excitation gap 
being larger or even one order of magnitude lager than 
the proximity-induced superconducting gap. The unique 
feature is in strong contrast to the case of the effective 
p-wave superconductors where the bulk excitation gap 




FIG. 4: (Color online) Numerical results of bulk excitation 
gap £gap of the TSC from the QAH state as a function of the 
proximity-induced superconducting gap A at /x = 0 (a) under 
different A with R, = 6 meV and (b) under different Vz with 
A = 4 meV. 


is always smaller than the proximity-induced supercon¬ 
ducting gap. This also applies to other QAH systems 
as seed materials to obtain robust TSCs and Majorana 
modes. 

Finally, we address the experimental feasibility of 
the TSC from the QAH state. Single layer graphene 
on the (III) surface of an antiferromagnetic insulator 
BiFe 03 can have an exchange field {Vz = 142 meV) and 
Rashba SOC (A = 1.4 meV), realizing a QAH insulator 
with a gap being 2Eq = 4.2 meVi^ This QAH state 
(IpI < Eq) in proximity to a conventional s-wave su¬ 
perconductor (eg., Nb with a large superconducting gap 
Anu = 0.83 meV— ) becomes a TSC since the topologi¬ 
cally nontrivial condition A^ + fj? < is easily satisfied 
due to |/x| < £0 <C |Ez| and |A| < IAnuI < |K|- With 
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A = 0.5 meV (A 2 = 0.21 meV < A < Ai = 2.1 meV) 
for estimation, we have the bulk excitation gap i?gap = 
2.05 meV, 1.56 meV and 1.06 meV, corresponding to a 
temperature of 23.8 K, 18.1 K and 12.3 K, at ^ = 0, 
0.5 meV and 1 meV, respectively. The large excitation 
gap (of the order of 10 K) ensures that robust Majorana 
modes can be achieved. 
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Appendix A: 14e(k) in Eq. (3) near the Dirac points 

Near the Dirac points, i.e., K = (47r/3,0) (r = 1) 
and K' = (—47r/3,0) (r = —1), iJe(k) in Eq. (3) can be 
expanded as a low energy effective Hamiltonian 


i7f(k) 


/ 14 

Vfirkx + iky) 

0 

\ iA_R(T-l) 


Vfirk^ - iky) 

14 

i\R(l + t) 

0 


0 

-*Afl(l -f r) 

-14 

Vfirkx + iky) 


*4_r(1 — 
0 

Vfirkx - 

-14 



(Al) 


The energy spectrum of this effective Hamiltonian 
is shown in Fig. 1(b). The minimum (maximum) 
energy of the conduction (valence) band is Eq = 
|14Ak|/ \/V^ + {—Eq) after a simple calculation and 

then the band gap is given by 2Eq. 


Appendix B: Gap closing condition of the BdG 
Hamiltonian HBdG(k) 

The gap of H^doi^) closes at the momenta P (single 
one), M (three inequivalent ones) and K (two inequiva¬ 


lent ones) points. Specifically, at the momentum P, the 
Rashba SOC vanishes [see Eq. O] , which is similar to the 
previous studies in semiconductors The gap closing 
condition is given by (/r±3t)^ = 14 ^ —A^ with -|- (—) rep¬ 
resenting lower (higher) energy band at P after a simple 
calculation. As for the momentum M, the Rashba SOC 
does not cause spin splitting but lead to an energy shift 
for the spin degenerate bands. We take M = (0, ^^^) 
for example and HBdG(Af) [see Eq. (3)] reads 



/ + 14 

%/3i-l + 

2 

0 

— A (2 + V^) 

A 

0 

0 

0 \ 



2 ^ 

—fj, -f 14 

-X{V 3 -i) 

0 

0 

A 

0 

0 



0 

—A(z + a/S) 

- 14 

2 ^ 

0 

0 

A 

0 


HudciM) = 

-X{V 3 -i) 

A 

0 

0 

2 '' 

0 

-M -14 

0 

0 

M + 14 

1 

1 

0 

0 

A 

A(? + y/S) 

(Bl) 


0 

A 

0 

0 

■y/Si+l -A 

2 ^ 

1^ + 14 

X{V 3 -i) 

0 



0 

0 

A 

0 

0 

X{i + a/S) 

1^-14 

2 ^ 



V 0 

0 

0 

A 

XiV^-i) 

0 

-A 

2 ^ 

/ 
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Performing a unitary transformation as EsdciM) = Uj^HBdG{M)UM with 


/ 


Um = 


V2 


0 

0 

2 

1 

0 

0 

0 

0 


one obtains 


HBdG{M) = 


f -t- fi-Vz 
0 
0 

2iX 

A 

0 

0 

V 0 


l — VSi 

2 

1 

0 

0 

0 

0 

0 

0 


0 

-t - -f 
2iX 
0 
0 
A 
0 
0 


14 


0 

0 

-1 + V3i 

2 

1 

0 

0 

0 

0 


0 

-2iX 

t- fi-Vz 
0 
0 
0 
A 
0 


-1+V3i 

2 

1 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

1-V3i 

2 

1 


0 

0 

0 

0 

l — y/3i 

2 

1 

0 

0 


-2iX 

0 

0 

t — fi Vz 
0 
0 
0 
A 


A 

0 

0 

0 

t ^ — Vz 
0 
0 

-2iX 


0 

0 

0 

0 

0 

0 

-l+Vdi 

2 

1 


0 

0 

0 

0 

-1 + V3i 

2 

1 

0 

0 


0 

A 

0 

0 

0 

i M + 14 
-2iX 
0 


0 

0 

A 

0 

0 

2iX 

— t + ^—Vz 
0 


(B2) 
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At ^ t, the block with the diagonal terms being —t — /r =F 14 and t + /r =F 14 in HBdG{M) is far from gap closing 
whereas the gap closing is determined by the remaining one. By considering that |A| <C t, we use the Lowdin partition 
metho d'^^d^ to obtain the effective Hamiltonian for the block determining the gap closing as 




^ t - ^ 

V 


14 + ^2^ 0 ^ A 0 \ 

0 f — ^ +14 + 0 ^ 

A 0 ^ -t + ^,-Vz-^ 0 

0 A 0 —t + fi + Vz + -^+Vz ' 


(B4) 


Then, the gap closing condition is (t — /r — 14 + M — 14 — — A^ = 0 or (f — ^ + 14 + M + 

14 — = 0- As |141 <C t, both conditions become {t — fj, + — A^ approximately. Furthermore, by 

considering that |A| <C t, we neglect the energy shift of 2X? jt and then the gap closing condition at the momentum M 
with /r ^ t is given by {t — /r)^ = — A^. Similarly, the gap closing condition at M with /x ~ —t is (t + /x)^ = — A^ 

under the approximation |A|, \Vz\ -C t. 

In contrast to the momenta F and M, the Rashba SOC at the Dirac points contributes to a finite spin splitting. 


Specifically, with K = (47r/3, 0), i7BdG(Ar) [see 

Eq. (3)] 

can be written as 
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which can be divided into two independent 4x4 parts, i.e.. Hi {H 2 ) without (with) the Rashba SOC terms. Then, 
we have 


/ +14 
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0 

- Vz 
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A 
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/X + 14 
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V 0 

A 

0 

H-Vz 


(B6) 


which is exactly the same as the Hamiltonian of semiconductors with the Rashba SOC, magnetic field and proximity- 
induced s-wave superconductivity at the momentum Fj^d^ This indicates that both have the same gap closing con¬ 
dition, i.e., = fj? + As for H 2 (not shown), due to the existence of the nonzero Rashba SOC terms, the 



























gap is always opened. Therefore, the gap closing condition at K is just the one in Hi part. Similar analysis can be 
applied to K' and we obtain the same gap closing condition as K. _ 


Appendix C: Numerical method for calculating 
Majorana edge states in zigzag and armchair 
graphene ribbons 

We investigate the Majorana edge states near the Dirac 
points in both zigzag and armchair graphene ribbons. We 


first study the case of zigzag configuration. The Hamil¬ 
tonian of zigzag ribbon can be obtained from Eq. (1) by 
choosing a unit cell and performing a Fourier transfor¬ 
mation along the direction parallel to the edge (assuming 
x-direction). Note that the unit cell of the zigzag ribbon 
is the same as the one in Ref. H^. Specifically, 


H 


zigzag 


[1 + \sgn{Xj^ - -f ^ ^ (ct14 - 

fex {h,32)iy fex JO' 

+ AE + H.C.) + zAE E 


fex 3 fex <jl,j2>OCT' 




(Cl) 


r 


where Xj^ — Xj-^ is the relative position between j 2 -th and 
ji-th atoms in the unit cell along the x-direction and 
sgn stands for the sign function. By exactly diagonaliz¬ 
ing i?zigzag) one obtains the eigenvalues and eigenstates. 
However, this method fails due to the computational limi¬ 
tations when the width of the ribbon becomes very large 
(eg., of the order of 10^ atoms in the unit cell in our 
calculation). Alternatively, the zigzag ribbon with the 
leading term, i.e., the hopping term, can be solved ana¬ 
lytically near the Dirac points4^ Near AT (r = 1) and K' 
(r = — 1), the eigenstates are given by 

_ ^gi('r|^S'H-fex)a: 

-Vf[{z - Tkx)e^y + iz + Tkx)e~''y]/e \ 

gzy _ g-zy I 

I 


with the eigenvalues being = v‘j{k'^ — z^) and A = 

\J | 2 (e^'^-^-e^-^)/z- 8 L| ' ^ width of the ribbon and 

z is determined by the equation = {kx + tz )/ {kx — 

Tz). Note that if zq is a solution of this equation, so 
does —zq. As only one of these two 

equivalent eigenstates needs to be taken. Then, one can 
use these eigenstates in Eq. (|C2I) with additional spin 
and particle-hole degrees of freedom included to con¬ 
struct complete basis functions for i/zigzag- We diago¬ 
nalize the Hamiltonian matrix of i?zigzag and obtain the 
energy spectrum and wavefunctions as shown in Fig. 2. 

We turn to the case of armchair graphene ribbon with 
the Hamiltonian being 


.^armchair — ^ E E xJO- + aEE( +H.C.) 

fcx <jlj 2 ><T fcx JCT fex 3 

+ *^E E '^h32~’^y ^.?1.72)4x7ig''^fcxj2<T' ~ ^E E ^^“’(^3* ,32 + '^jT-H,Jz)4x7i*g'‘^fcxJzg' + H.C.] 


fcx (31,32)<T<t' 


+ *AE E + |<"')]cU<zCfexj5^' +H.C.} (C3) 


Ji 32 


kx 3 i32^ 

2 


in which j* (j^) represents the j J-th (j|-th) atom of the armchair ribbon with only the hopping term analytically 

first (fourth) column in the unit cell. Note that the unit 

cell of the armchair ribbon is the same as the one in 

Ref. 113 and the edges lie along the x-direction. Similar to 

the case of the zigzag graphene ribbon, we first solve the 
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near the Dirac points. The eigenstates read 


= 2Ae*'=^-sin[(|i^| + fc„)y] 

^ ^ -Vfikoo - ikn)/e ^ 


(C4) 



y 

FIG. 5: (Color online) (a) Energy spectrum of armchair 
graphene ribbon in the presence of the Rashba SOC, exchange 
field and s-wave superconductivity from the proximity effect, 
(b) ((c)) Real space probability amplitude \ip\ across the width 
for the Majorana edge state with a smaller (larger) momen¬ 
tum \kx\ {kx < 0) at one edge (i.e., y = 40000) (only part 
of the ribbon is shown). The fluctuations of li/il at the posi¬ 
tions far away from the edge are due to numerical error. Here, 
14 = 6 meV, A = 4 meV, y, = 2 meV and A = 2 meV. 


where the eigenvalues are = v‘j{k'^ + k"^) with kn = 
mrlL — \K\ and A = These eigenstates construct 

complete basis functions for i?armchair with additional 
spin and particle-hole degrees of freedom. By diagonal¬ 
izing the Hamiltonian matrix, one obtains the energy 
spectrum and eigenstates of armchair graphene ribbon 
as shown in Fig. [5] In Fig. EKa), we find that there exist 
eight zero energy states, corresponding to four Majorana 
fermions at each edge, which is similar to the case of 
zigzag ribbon. We then show the real space probability 
amplitude of two Majorana edge states at the same edge 
with a smaller and larger momentum \kx\ {kx < 0) in 
Figs.EDb) and (c), respectively. It is seen that both show 
obvious decays and oscillations but the decay lengths and 
oscillation periods are different. 
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